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Abstract 

The projection pursuit index defined by a sum of squares of the third and the fourth 
sample cumulants is known as the moment index proposed by Jones and Sibson fl4|. Lim- 
iting distribution of the maximum of the moment index under the null hypothesis that the 
population is multivariate normal is shown to be the maximum of a Gaussian random field 
with a finite Karhunen-Loeve expansion. An approximate formula for tail probability of the 
maximum, which corresponds to the p- value, is given by virtue of the tube method through 
determining Weyl's invariants of all degrees and the critical radius of the index manifold of 
the Gaussian random field. 

Key words: critical radius, Euler characteristic heuristic, Hotelling-Weyl tube formula, max- 
imum of a Gaussian random field, multiple testing, sample cumulant. 
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1 Introduction 

1.1 Assessing the significance in projection pursuit 

Suppose that for each of n individuals, a q dimensional random vector G M'^, t = 1, . . . , n, is 
observed as an i.i.d. sample. In the analysis of such multidimensional data, the projection of the 
q dimensional data onto a lower dimensional subspace is often used for the sake of interpreting 
the data. In such cases it is important to select the subspace which clarifies features of the 
data interesting to the analyst. In the principal components analysis or the canonical correla- 
tion analysis, the subspaces are selected based on the variance of data ([5]). The exploratory 
projection pursuit is the method for detecting the subspace based on the non-normality of data 
(|12j). As a similar method, the Fast ICA (independent component analysis) is known (yjj). 

Let §'^~^ be a set of q dimensional unit vectors, or the set of directional vectors in W^. In the 
one dimensional projection pursuit, for each directional vector h G S'^^^, the projection pursuit 
index /ri(/i) is defined as a measure for the non-normality of one dimensional projected data 

zt = {xt,h) eM., t = l,...,n, (1.1) 

and then the direction h* = argmax/„(/i) attaining the maximum of the index is searched. 

However, the index Inih) is a random function of h depending on the samples xt^s. Even 
when xt's were distributed according to the multidimensional normal distribution, the function 
Inih) is not constant, and the direction h* which achieves the maximum exists. Therefore, it is 
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important to assess whether it is not a pseudo peak caused by stochastic fluctuations. For this 
purpose, the framework of the multiple testing can be employed. Consider the null hypothesis 
that the data are distributed according to the multidimensional normal distribution 

Ho : xt ~iVg(/i,S) i.i.d., (1.2) 

and let 

Fn{c) = P( max >c\Ho) 

be the upper probability of the maximum of /n(/i) under the null hypothesis. Then, 

is the p-value in the sense of multiple testing, and we can use the p-value as a measure of the 

significance of the maximum (Sun |2Uj). 

Sun [20] described the limiting null distribution of the maximum of Friedman [8]'s projection 
pursuit index in terms of a Gaussian random field as sample size goes to infinity, and gave an 
approximation formula for it by an integral-geometric method referred to as the tube method 
( ) • In this paper we gives an approximation formula for the moment index proposed by Jones 
and Sibson [H] by the tube method. 

The moment index treated here is as follows: Let Kk^n{h) be the kih sample cumulant of the 
projected data zt in (fTT]) . and let Bi^n{h) = K^,n{h) / K2,n{hf^'^ and B2^n{h) = i^4,n(/i)/-f^2,n(/i)^ 
be the sample skewness and the sample kurtosis, respectively. Then the moment index is defined 
by 

In{h) = ^B^,^{hf + ^B2,n(.h?. (1.3) 

Differently from Sun [2D]'s treatment for Friedman's index, we can determine geometric invari- 
ants of all degrees, and accordingly give an accurate formula for the p-values. 

The structure of the paper is as follows: The main results are summarized in Section [2l 
There, the limiting distribution of the maximum of the moment index is described in terms of 
a Gaussian random field with a finite Karhunen-Loeve expansion, and determine the geometric 
invariants of the index manifold. An approximation formula for the upper probability of the 
maximum can be obtained by incorporating these invariants. Some numerical experiments to 
examine their accuracy are given there. The main results of Section [2] are proved in Section [H 
Prior to Section HJ we give a brief summary of the tube method in Section [3] as far as required. 

1.2 The tube method 

Here we give a very brief historical review of the tube method. 

As explained in Section [3l the term tube means a spherically tubular neighborhood around 
a set in the sphere. Hotelling [10] pointed out a relation between the p-value of a testing 
problem in nonlinear regression and the volume of the tube, and demonstrated to calculate the 
p-value by presenting a one dimensional formula for the volume of tube. Weyl |25j generalized 
Hotelling [IHl's formula to the general dimensional case. More recently, Knowles and Siegmund 
|15j and Sun [21] found out the relation between the formula for the volume of tube and the 
tail probability formula for the maximum of a Gaussian random field. Since then, the tube 
method was applied to statistical problems such as calculating null distributions of max-type 
test statistics, or adjusting the multiplicity in multiple testing problems. For example, the 
asymptotic distribution of the Anderson-Stephens statistic ([6]) for testing the uniformity of 
direction can be evaluated ([IS]). For the other examples, see [IT] and [TBJ. Nowadays, the tube 
method was proved by Takemura and Kuriki [22J to be a special case of the Euler characteristic 
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heuristics, which is known as another approach for approximating the distribution of the maxima 
of random fields developed by Adler ([2], [3]), Worsley (|26j, \27j) and Taylor ([23]). For recent 
developments of the tube method and the Euler characteristic heuristic, see Adler and Taylor 
[1]. See also [E]. 

2 Main results 

We begin with giving the limiting distribution of the moment index /n(/i) in ()1.3p under the 
null hypothesis of multivariate normality. Without loss of generality, we assume that xt^s are 
distributed according to the q dimensional standard normal distribution Nq{0,lq). 

Theorem 2.1 Let .^i G M^^, ^2 £ I^^* be random vectors consisting of independent standard 
normal random variables. For a unit vector h £ S'^^^, let 

Zi{h) = {h0h^h,^i), Z2{h) = {h0h0h0h,^2), (2.1) 

where (gi denotes the Kronecker product. Under the null hypothesis Hq in p.2|) . as n — > 00, 
max^g§q-i Inih) converges in distribution to max/^^gg-i I{h), where 

I{h) = Ziihf + Z2ihf. 

Proof Let C(S'^~^) be the Banach space of real valued continuous functions on S'^~^ endowed 
with the supremum norm. Note that the sample cumulant iCfc^„(-), the sample skewness i?i^„(-), 
the kurtosis -B2,n(")> ^^'^ the moment index /n(") are the elements of C(S'^~^). Theorem 2.1 of 
Kuriki and Takemura [17J states that the (■y/ni?i^„(-), \/n-^2,n(')) converges to {\^.Zi{-), \/4l-^2(')) 
in distribution in the space C(S'^~^). The theorem follows from the continuous mapping theorem. 



By means of Theorem 12.11 above, for large sample size n, the p- value Fn{In{h*)) can be 
approximated by F{In{h*)) with 



F(c) = P( max 1(h) > cV 



Moreover, by letting 



Z{h,e) = cos9Zi{h) + sm9Z2ih), h £ Be (-|,| 



with Zi{h) and Z2{h) given in (j2.ip . we have 

\ max 1(h)} = max Z(h,e). (2.2) 

UeS9-i J ?teS9-i.6»ef-7r/2.7r/2l 



Therefore, from now on, we restrict our attention to the distribution of the maximum of Z{h, 9). 
Let 

p = q^ 

and 

M=^{cos9{h(^h(^h),sm9{h(^h^h(g>h))eW\heS''-^, }. (2.3) 
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The maximum (|2.2p can be rewritten as 

max(x,0, e = (ei,---,ep)~A^p(0,/p). (2.4) 

Note that M and S"?"^ x (— 7r/2, 7r/2] are one-to-one. It is easy to see that M is a g dimensional 
closed submanifold of S^^^. As shall be explained in Section [Sj (12. 4p is of the canonical form of 
the tube method in (j3.ip . 

The upper probability function of the chi-square distribution with i/ degrees of freedom is 
denoted by 

The volume of the m — 1 dimensional volume of the unit sphere W^~^ is denoted by 

The following is the main theorem of this paper. The proof is given in Sections 14.11 and 14.21 
Theorem 2.2 As c ^ oo, 

e=0, e:even 



where 



^ _^ (-3)-/^(g-l)! ^ (g-e-2j) 



and 



/.7r/2 

£;fc= / (3cos2e + 4sin2 0)'=d0 (2.^ 



Remark 2.1 ii^jfc in p.Sp with A; an integer or a half- integer can be evaluated numerically by 
recurrence formulas: 

7(2A;-1) 12(A;-1) , , 3 , , 

Ek = ^^, ^^fc-i S for A: = l,-,2,..., 2.10 

Ik k I 

and 

7(2A; + 3) fc + 2 ^ = _ ^ 2, . . . , (2.11) 

24(/c + l) 12(/c-M) 2' ' 2' ' ^ 

with the boundary conditions 

^^/2 = 4£;(l/4) = 4 X 1.46746, ^0 = tt, ^_i/2 = i^(l/4) = 1.68575, E^x 



vr 



2^/3' 

where £'(1/4) and if (1/4) are complete elliptic integrals of the first kind and the second kind 
([1], p. 608-9). The proofs for (f2Tn]) and dTTTI) are given in Section Ol 
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To conclude this section, we give numerical examples for the purpose of examining the 
accuracy of the formula. The tail probability of the maximum for g = 2 is given by 

P(^max /(/i) > c^) ~ w{G3{c^) - Gi{c^)} = w^ce'''^''^, c ^ oo, (2.12) 

where w = 2E{l/A) = 2 x 1.46746. 

Figured] depicts the empirical upper probability of the limiting distribution P(max/igg2-i I{h) 
> x) estimated by Monte Carlo simulations based on 10,000 replications, and its approximation 
by the tube method. One can see that the quantiles of the limiting distribution are fully 
approximated by the tube method approximation (j2.12p . 

Figure [2] depicts the empirical upper probability of the finite sample distributions P(max;^gg2-i 
In{h) > x) when n = 300, 1000, 3000, oo. The number of replications is 10,000. 




- tube approx. 



Figure 1: Tail probability of limiting distribution (solid line) and its approximation by the tube 
method (dotted line). 



3 Summary of the tube method 

3.1 Volume of the tubes and tail probabilities of the maxima 

In this section we summarize the facts on the tube method required for proving Theorem 12.21 
We state Theorem 13.11 since its statement is not given in existing literature. 

Let be the unit sphere in W, and let M be a closed subset of S^~^. Assume that M is 
a d dimensional closed submanifold without boundaries embedded in S^"^, and is endowed 
with the metric induced by the standard inner product (•, •) of W. 

The set of points of S^"^ whose great circle distance (angle) from M is less than or equal to 
a constant 9 is called the tube about M with the radius 6, and denoted by 

Tube(M,6') = \y e S^"^ | dist(y,M) < o}, dist(y,M) = min cos~^(y,x). 

In a similar manner, the Euclidean tube is defined in the Euclidean space by the usual distance. 
But it does not play any role in this paper. 
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Figure 2: Tail probabilities of finite sample distributions (n = 300, 1000, 3000, oo). 



Let y be a point of S^^^M. The point x = pr(y) which attains the minimum min^^g m dist(?/, x) 
is called the projection of y onto M. If y is close to M, then pr(y) exists uniquely. Whereas, if 
y is far from M, then there can exist two points xi,X2 G M equidistant from y which attain the 
minimum min^j^M dist(?/, x) simultaneously. The supremum of the distances which assures the 
uniqueness is called the critical radius. 

Definition 3.1 When the pr(y) G M is defined uniquely for every y S Tube(-M, ^) \ M, it is 
said that the tube Tube(M, 9) does not have a self-overlap. The supremum 

6c = sup{9 > I Tube(M, 9) does not have a self-overlap} 

is called the critical radius of M. 

The volume of a tube whose radius is less than or equal to the critical radius 9c can be calculated 
by taking a coordinate system based on the projection (the Fermi coordinates). The following 
proposition for the dimension d = 1 is due to Hotelling [10] , and due to Weyl [25] for the general 
dimensional case. Here Q,p denotes the p — 1 dimensional volume of defined in ()2.6p . and 



T{a + b) 



(1 - t)^-^dt 



T{a)T{b) 

is the upper probability function of the beta distribution with parameters (a, b). 
Proposition 3.1 For < 9 < 9c, p — 1 dimensional volume of the tube is given by 

d 

Vol(Tube(M,^)) = 17p ^ i^eJe{0), 



e=0, e:even 



where 



e)/2(cOS^^), 



nn) r((d + i-e)/2) 

and the Ke is the intrinsic invariant of the manifold M defined below in (13. 6p . referred to as 
Weyl's curvature invariant. 
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Let ^ = (^1, . . . , ^p) be a random vector consisting of independent standard normal random 
variables. That is, ^ ~ Np{0,lp). Define a Gaussian random field on a submanifold M of EP~^ 
by 

Z{x) = {x,0, xeM {CSP^'^). (3.1) 

This is a canonical form of Gaussian random fields of mean and variance 1 with a finite 
Karhunen-Loeve expansion. 

By replacing QpB by the upper probability function of the chi-square distribution G in (j2.5p . 
we have an approximation formula for the tail probability of the maximum of Z{x) (|l7], [22j ) . 

Proposition 3.2 As c —>■ (yo, 

d 

l+tan2ec)c2/2) 

e=0, e:even 

where 



, , , r((d+l-e)/2) - . 2, 
Note that the larger the critical radius 6c is, the smaller the order of the remainder term is. 
3.2 Weyl's curvature invariants 

As we saw in Propositions 13 . 1 1 and 13.21 the Weyl's curvature invariants Kg and the critical radius 
9c of the manifold M are needed in applying the tube method. We will explain the way to 
determine them in this and subsequent subsections. 

Write a local coordinate system of a d dimensional closed manifold M as (t*). The metric 
tensor is denoted by gij, and write the (z,j)th elements of the inverse of the d x d matrix [gij) 
as g^K Abbreviate d/df to di. The connection coefficients and the curvature tensor are given 

by 

d ^ 
^ij = Yl 9''''^ij,h with r,j- fc = -{digjk + djgik - dkgij), (3.2) 

h=l 

and 

d 

R'lj = Y,g'^R\,j with R'k^,=^,T'^,-^JTi, + Y,{^lT%-T'^,^|,), (3.3) 

h s=l 

respectively. Let 

H^ = Rl^j- {5^5^-616';), (3.4) 
where is Kronecker's delta. For e = 0, 2, . . . , [d/2] x 2, let 

Here the summation Y2i taken over all sets of e/2 paring made of distinct elements of 
{l,2,...,d}, that is, all possible ways of {ii,i2, ■ ■ ■ ,ie} C {l,2,...,d} satisfying ii < 12, 
is < i4,...,ie_i < ie and ii < is < ■■■ < ie-i- The summation Y2ct taken over all per- 
mutations a of {1,2, ...,e} such that a{l) < cr(2), a{3) < (t(4), . . . , cr(e — 1) < cr(e). Then, 
Weyl's curvature invariants are defined by 

Kc= I He det{gij)^/^dt^ ■■■dt'^, e = 0,2, . . . ,[d/2] x d (3.6) 
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(Weyl [25]). 

For instance, for e = 0, 2, 4 are given as follows: Hq = 1, and hence kq is the d dimensional 
volume of M. 

^= E 4 = ^E4 = ^{E4-'^('^-i)}' 

l<i<j<d i,j=^ iii=l 

where Ylij=i scalar curvature. 

l<i<j<k<l<d 

Tjij Tjkl I Tjjfc Tjjl Tjil Tjjk Tjjk TTll , Tjjl Tjlk Tjkl Tjij 



i j',fc,/=l 
^ a! 2 <^ 

8{(E-^ii) ~^ E E ^ij-^M 

i,j=l i,j,k,l=l i,j,k,l=l 

d 

-2{d - 2){d - 3) E + d{d - l){d - 2){d - 3)}. 



See Gray ([9], Lemma 4.2) for the invariants of a Euclidean tube. 
3.3 Evaluation of critical radius 

In this subsection we give theorems useful in calculating the critical radius of a closed subman- 
ifold of the sphere. 

Proposition 3.3 The critical radius 9^, of a closed submanifold M of satisfies 

cot^ 6'c= sup h{x,y), h{x,y) = ^ — T^AtI' (3-7) 

where Px is the orthogonal projection onto the linear subspace spanjx} © T^M of M^, and TxM 
is the tangent space of M at x ( [13] , [IZj ) . 

A theorem corresponding to a Euclidean tube is given by Federer ([J|, Theorem 4.18). 
The radius 6^°'^ satisfying 

cot'^0^°''= limsup h{x,y) (3.8) 

is called the local critical radius, which is characterized as the curvature radius of M at x (|13j. 
|17j). By definitions, 9^°'^ > 9c, and the equality holds if the supremum in ()3.7p is attained when 
||y — — > 0. 

Define a real-valued function on M x M by r(x,y) = {x,y). This is the covariance function 
of the Gaussian random field (j3.1|) . Denote the local coordinate system about x and y by (s*), 
(f), respectively. 
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The set of the critical points of r{x, y) which are not contained in the diagonal set is denoted 

by 

C = e M X M \ X ^ y, —r{x,y) = 0, —j{x,y) = o|. 

Then we have the following theorem. 



Theorem 3.1 The critical radius 9c satisfies 

9c = niin 



{9'r inf icos-i(x,y)}. 

(x,y)&C I J 



Proof By Lemma 5.2 of [23], if the supremum of h{x,y) is attained at a point not contained 
in the diagonal set, then it belongs to C. Furthermore, for the points {x, y) £ C, it holds that 
PxV = {x,y)x, 

,1 N l + {x,y) ^(^ -i, C\ 

and hence ^ 

sup y) = cot^ ( inf - cos^"'^(x, y) ) . 
{x,y)ec ^{x,y)<^C 2 / 

Since the supremum of h{x, y) over the diagonal set is cot^ the theorem follows from Propo- 
sition 13. 3[ I 



A theorem corresponding to a Euclidean tube with the dimension c? = 1 is given by Johansen 
and Johnstone ([13], Proposition 4.2). 



4 Proof of Theorem 12.2 
4.1 Proof of (Em 

In this section, we prove Theorem 12.21 By means of Proposition 13.21 the approximation formula 
for the upper probability of the maximum can be given through determining Weyl's curvature 
invariants Kc and the critical radius 9c of the index manifold M in (12. 3p . The former is given 
here, and the latter is given in the next subsection. 

The metric tensor, the connection coefficients, and the curvature tensor for M are denoted 
by g, r, and R, respectively, as in Section 13.21 Also, the same quantities for S''"^ are denoted 
by g, r, and respectively. 

Write an element h of §^~^ by a local coordinate system as h = h{t), t = (f). Let hi = dh/df. 
The metric of S*^"^ is gij = {hi, hj). 

An element x of M can be written as 

X = {cos 9{h (g) /i (g) h),sm9{h (g) /i (g) /i O /i) ) € M 
bases of the tangent space of M are 

= (cos 9{hi ®h®h + h®hi®h + h®h® hi), 

sin 9{hi ®h®h®h + h®hi®h®h-\-h®h®hi®h 
+h (g) h (g) h ^ hi)), i = l,...,q—l, 

= {- sin 9{h®h® /i), cos 9{h®h®h®h® h)). 



in terms of {t,9). The 
dx 



dx 
89 
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In the following, 9 is regarded as the 0th coordinate of t. The metric tensor of M is 



9i3 



v{e)g^,{t) ifi,i/0, 
1 if i = j = 0, 

otherwise, 



where 

v{e) = 'icos^ 6 + 6 = 2, + 9 = A - cos^ e. (4.1) 
From this, the volume element of M is shown to be 

Aei{gij{t)Y/^dt^ ■ ■ ■ dt"'^ (4 2) 

Note that det{gij{t))'^/'^dt^ ■ ■ ■ dt^-^ is the volume element of S^"^ 

Let V and v be the first and second derivatives of v = v{9). After some calculations along 
the lines with (13.21). it is shown that the non-zero connection coefficients of M are 



r^, = f f^, rfo = rgi = ^^5^ r". = -^m, (ij, k + o). 



and all of the other coefficients are 0. 

Next we will derive the curvature tensor by ()3.3p . Put 



Jij = Oi Oj - didj . 

Noting that the curvature tensor of the unit sphere S^^^ is R^j = Jf-, after cumbersome calcu- 
lations we see that the non-zero elements are 



R 



.V 4\v 



Furthermore, noting that v = 2 cos sin 0, (w)^ = 4 cos^ 9 sin^ 9 = 4(4 — v){v — i) = — 4(f ^ — 7v + 
12), V = 2 cos^ 9 — 2 s\v? 9 = 2(4 — v) — 2{v — 3) = —2(2?; — 7), we have the non-zero elements of 



H^l in ^ as 

H^j = a4J, = -H^^ = -Hl^ = i/Of = /3<5f (.,i, A;,/ / 0), 

where 

6 12 12 

a = a{9) = — + ^, (5 = (3{9) = — ^. 

V V V 

We substitute these quantities into (|3.5|) to obtain He, e = 0, 2, . . . , [g/2] x 2. 

(i) The case where the set of the indices ii,i2, ■ ■ ■ ,ie in the right-hand side of (13. 5p does not 
contain 0. Because the number of the ways to make e/2 pairs from q — 1 distinct objects is 



{q-iy. 



(g-l-e)!2^/2(e/2)!' 
the summation of all terms corresponding to the case (i) becomes 

" ^,-l-e)!2e/2(e/2)!- ^'•'^ 

(ii) The case where the set of the indices ii, • • • , in the right-hand side of ()3.5p contains 
0. In this case, ii = 0, and i2 / 0, a{l) = 1 (v(i) = 0). Noting that there are q — I ways for 
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12 (^2 = 1, • • • 1 — 1)) and that is, ^4, . . . , are indices resulting from making e/2 — 1 pairs from 
the set {1, 2, . . . , g — 1} \ {^2} having q — 2 elements, the summation of all terms corresponding 
to the case (ii) becomes 



Summing up (14. 3p and (14. 4p along with 



a 



and 



yields 



where 



'JX A- 



j=0 



(-3r/2(g-l)! 

(g-e-l)!(e/2)!' 



and for j / 0, 



' ^ ^ ^ U(9-l-e)!2^/2(e/2)!V i ; (g - e)!2-/2-i(e/2 - 1)! V i - 1 7 

(_3)e/2(g_i),(g_g^2j)(-2y 

(g-e)! {e/2-j)\j\ 
Since the expression for Aj with j 7^ is consistent with that for Aq, we have 

_ (-3r/2(g-i)! ^ (g-e-2j)(-2y 1 

(g-e)! (e/2-j)!j! W^+r 

Finally we obtain Kg in (|2.7|) by integrating Hf, over M with respect to the volume element (|4.2p . 

4.2 Proof of dan]) 

In this subsection, making use of Theorem 13. H we show that the critical radius of the index 
manifold M in ([231) is Oc = tan-^(3/4). This implies that pc = l + tan^ 6^ = 25/16. Throughout 
this subsection, we assume that vectors are column vectors for notational convenience. For 
instance, (x, y) = x'y, where ' denotes the transpose. 

We begin with obtaining the local critical radius 9^°'^ by ()3.8p . Let 

cos 9{h h) \ cos6{h®h®h) 

sm.9{h®h®h®h))' ~ \sm.9{h®h®h®h) 
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be two points of M. Write for simplicity hi = dh/df, Xi = dx/df, xq = dx/dO, and G = (gij), 
V = v{9) defined in ()4.ip . The orthogonal projection matrix onto spanjx} © TxM is denoted by 
Px- Since span{x} is orthogonal to T^M, we have 



X -P^ — ^>-C X ^ X X • • • ^ X X q — l;*-^ "^o) 



vG 



(x x) + {x Xi, . . . ,X Xq^i){vGy 



/ x'x\ 

\x'Xq-X, 



I x'x \ 
x'xi 

x'Xq^l 
\ x'xq J 

+ [x'xof 



The first term of the right-hand side is the square of 

x'x = {h'hf cos ecose+ {h'hf sin OsmO. 

Noting that 

x'xj = wh'hi, w = 3(/i'/i)^ cos ^ cos + 4(/i'/i)^ sin 6^ sin 0, 
the second term becomes 

/ h[ 



— h'{hi,... ,hq-i)G 



-1 



w 



h' = — h'ilq - hh')h = — (1 - (h'hf^ 

V V 



The third term is the square of 

x'xo = - {h'hf cos 61 sin 6* + (h'h)^ sin fl' cos 6*. 
Summing up these three terms, the numerator of the right-hand side of (j3.8p is 



w 



1-xPxX = 1 - {{h' hy COS e COS 9 + {h' hf sin 9 sin ey- — {l-{h' h) 

u\3 „„„ 2 a ( fu' u\4: „■„ Q „„„ o\2 



{-{h'hy COS 9 sin 9 + {h'hy sin 9 cos 
- cos^ ^ cos'^ 9 — cos^ ip sin^ 9 
(3 cos'^ ijj cos 9 cos 9 + 4: cos^ ip sin 9 sin 9)"^ 



3 cos2 9 + 4 sin^ 9 



sin^ 



= / (say), 

where h'h = cosip. On the other hand, the denominator of the right-hand side of ()3.8p is 

(1 — x'x)^ = (1 — cos'^'i/'cos(9cos6' — cos^ 'i/'sin(9sin6')^ 
= g (say). 

The local critical radius 9^°^^ can be obtained by cot^ 9^°^ = limsup//^ when 9 — 9^0, 
■0 — > 0. Let 9 — 9 = 6 and u = sin^ 9. Ignoring i/;^, t/^^J^, and 6^ as infinitesimals, we have with 
aid of symbolic calculation that 



f^3il + u)i;^ + -^i;'5' and g ^ 

6 + u 4 
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Letting 5"^ ~ kifp' for a constant k (may be or oo), we have 

/ ^^ (l + ^)(3 + u) + 4fc 
5 (3 + n)(A: + 3 + n)2" 

As a function of /c, the right-hand side of the above takes its maximum 

48 

(3 + n)2(3-n) 

at /c = (3 + u){l — u)/2. Furthermore as a function of n, this takes its maximum 48/27 = 16/9 
at n = over < u < 1. Note that when u = 0, A; = 3/2 and = 0. 

Summarizing the above arguments, one can see that hmsup//^ = 16/9 is attained when 
61, 6* ^ 0, ^ 0, \0 — 6\ ^ y^3/2^/;, and accordingly 



e^^^ = cot-i(yi679) = tan-i(3/4) = 0.205^. 

As the second step, we confirm that the local critical radius is really the critical radius. The 
covariance function of (12. 4p is 

x'x = cos6 cos 6 (h'h)^ + sin 9 sin 6 (h'h)'^ 
= cos 6 cos 9 cos^ ip + sin 6 sin 9 cos^ V 
= r 0^,9, 9) (say). 

The ranges of the variables are 

VG[0,7r], 9,9 (4.5) 

The set of the critical points are the set of the solutions of 

~ ~ 

= — - = — sin-(/'(3 cos cos 6* cos^ V + 4sin sin0 cos^ V)) (4-6) 
dtp 

Qt' ~ 

= — = — sin^cos^cos^ V + cos6'sin0cos^ V) (4-7) 
o9 

= ^ = — cos sin 6^ cos^ '(/' + sin0 cos ^cos^ V- (4-8) 
d9 

(i) The case sin ^ 7^ 0. From (14. 6p . (3 cos 9 cos^ ip, 4 sin 9 cos^ V') is orthogonal to (cos 0, sin 0). 
From (j4.7p . (cos cos'^ "0? sin cos"^ ^) is orthogonal to (— sin 0, cos 0). Combining these, 

= 3 cos^ 9 cos^ -0 + 4 sin^ ^ cos'^ -0 = cos^ V(3 cos^ ^ + 4 sin^ 6^ cos^ 0) , 

from which cos0 = 0. Because of (14. Sh . = 7r/2. Conversely, when -0 = '7r/2, ()4.6|) - ()4.8p are 
satisfied. Thus, 

r = 0, - cos-i = ^ > 6*^^°^ 
2 4 

(ii) The case sin0 = 0. Then cos0 = ±1 (^ = 0, vr). In this case both (j4.7p and (j4.8p are 
reduced to sin(6' =F ^i) = 0. Because of gSD, f = If -0 = and 6* = 6i, then (/i, 6*) = (/i, ~9), or 
a; = X. Hence, it should be = vr, = —0, and 

r = -l, lcos-i(-l) = |>0l-. 

We have proved that the critical radius is attained locally. 
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4.3 Proof of the recurrences (EHOl) and 

For V = v{9) = Scos^ 6 + 4:sm^ 6 = A- cos^ 6 = 2, + sin^ 6*, 

'■Tv/2 



Ek = / v{efde 

/7r/2 />7r/2 
{A-cos^e)v^-^de = AEk-i- / cos^V^^d^ 
-7r/2 y-7r/2 

l-Tv/2 

+ / sme{cose?;'=-^}'d0 

y-7r/2 



'-7r/2 J-7r/2 

7r/2 ^7r/2 

4£'fc_i - sin 6* cos 6'?;''"^ 



-7r/2 J-7r/2 
f.7r/2 i-Tv/2 



and hence 



or 



sin^ ^v'^^^de + / sin cos e{k - l)v^-^2 sin cos OdO 

-tt/2 J-n/2 

/n/2 f-TT/2 
{v - 3)v''~^de + 2{k -I) I (t; - 3)(4 - v)v^-^de 
-7r/2 J-tt/2 

AEk-i -Ek + 3Ek-i - 2{k - l)Ek + Uik - l)Ek-i - 24(A: - l)Ek^2 
(_2A; + l)Ek + (14fc - l)Ek-i ~ 2iik - l)Ek^2, 

2kEk = 7{2k - l)Ek-i - 24{k - l)Ek-2 
-2{k + 2)Ek+2 + 7{2k + 3)Ek+i = 24{k + l)Ek. 
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